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i. 


ESTIMATION  OF  PARAMETERS  FOR  TRUNCATED  DISTRIBUTIONS 

1 .  INTRODUCTION 

The  analysis  of  chemical  warfare  systems  involves  the  use  of 
distributions  of  droplet,  particle,  spot,  and  fragment  sizes.  Many  field  and 
laboratory  tests  are  conducted  to  obtain  data  for  the  determination  of  such 
distributions.  As  a  practical  matter,  the  experimenter  usually  assumes  the 
data  are  best  characterized  by  a  lognormal  distribution,  which  is  equivalent 
to  assuming  that  the  logarithm  of  the  sizes  is  normally  distributed.  Then  he 
proceeds  to  perform  the  routine  calculations  for  the  mean  and  standard 
deviation  of  the  data  and  reports  the  results.  If  his  experimental 
measurements  span  the  range  of  possible  sizes  and  the  number  of  observations 
is  sufficiently  large,  this  procedure  leads  to  acceptable  results.  However, 
many  measuring  devices  are  limited  to  a  specific  operating  range,  particularly 
on  the  low  end,  beyond  which  no  data  are  obtainable.  Since  most  experimenters 
are  familiar  with  procedures  for  dealing  with  the  effects  of  sample  size  on 
their  estimates,  this  topic  is  not  addressed  further.  On  the  other  hand,  many 
experimenters  fail  to  recognize  or  consider  properly  the  effect  of  such 
missing  data.  The  purpose  of  this  report  is  to  describe  the  problem  and 
present  a  method  for  determining  the  best  estimates  of  the  parameters  of  a 
truncated  lognormal  distribution. 

2.  THE  PROBLEM 

For  the  present  purpose,  it  is  assuned  that  a  lognormal  distribution 
of  sizes  of  spherical  liquid  drops  is  adequate  to  represent  the  data  collected 
in  an  experiment  which  counted  the  number  of  drops  in  each  of  a  finite  set  of 
size  intervals.  The  mathematical  expression  for  the  probability  density 
function  of  the  lognormal  distribution  is 


f(d)  »  (din  a  /2ir)_l  exp 


1-1/2  iln,d:*n  dN. 


in  a, 


g 


where 


d  *  diameter  of  a  drop 

d^  »  number  median  diameter  (geometric  mean) 
o  -  geometric  standard  deviation  of  the  diameters 
g 

The  objective  of  the  experimenter  is  to  obtain  the  best  estimates  of  dN  and 
Og,  or  possibly  the  mass  median  diameter  d^,  which  is  related  to  d^  by  the 
expression 


in  dM  -  in  dN  +  3  in2 


(2) 
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Suppose  the  experimental  data  obtained  consist  of  the  number  n^  of 

drops  in  the  i-th  size  interval  such  that  d^_^  <_  d  <  d^  for  i«l,2,...,k  and 

d  >  0.  No  information  is  available  for  any  d  <  d  .  If  this  truncation  of 

o  o 

the  distribution  at  dQ  is  not  considered,  the  results  reported  will  be 


a  /> 

where  the  estimates  d  and  o  may  be  far  from  the  best  estimates.  Depending 

N  g 

on  the  degree  of  truncation,  the  estimate  of  d^  will  be  too  high,  while  that 

of  a  will  be  too  low  by  an  undetermined  amount, 
g 

3.  A  SOLUTION 

In  1908,  Pearson  and  Lee*  reported  an  algorithm  for  obtaining 
maximum  likelihood  estimates  of  the  mean  and  standard  deviation  of  a  singly 
truncated  normal  distribution.  Maximum  likelihood  estimates  nay  be  defined  as 
those  which  describe  the  distribution  of  a  random  variable  most  likely  to 
produce  the  set  of  observations  of  a  random  experiment.  In  many  instances, 
these  estimates  are  unbiased  and  of  minimum  variance;  but,  these  properties, 
while  most  desirable,  are  not  guaranteed. 

Tine  Pearson-Lee  method  consists  of  solving  the  equations  derived  by 
equating  the  first  and  second  sample  and  population  moments  about  the  point  of 
truncation.  If  the  point  of  truncation  is  not  known  precisely,  it  is 
recommended  that  a  point  slightly  less  than  the  lowest  observed  size  be 
used.  The  equations  to  be  solved  are 


*  Pearson,  Karl,  and  Lae,  Alice.  On  the  Generalized  Probable  Error  in 
Multiple  Normal  Correlation,  Blometrika  6,  59-68,  (1908). 
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and 


A  /A 

a  a 

-  - ' x| 
u  \u 


rk  2 

Zi-1  niUi 


,1  -  „! 


where  u 


original  random  variable 

truncation  point  (u1  <  u1  ) 
o  — 


*i-i  Vi 

Ek  n 
i»l  i 


I  A 

U1  -  a 


•,  a  standard  deviate 


where 


A 

ffi 


mean  of  the  normal  distribution  (estimated) 


*  estimated  normal  standard  distribution 


z(x)  «  (/2x)  *  exp  (-  x^/2) 


♦<x)  -  J1 


z(t)  dt 


(6) 

(7) 


(8) 


(9) 


(10) 

(ID 


By  solving  equation  (5)  for  a/u  and  substituting  in  equation^(6) ,  a  nonlinear 
equation  in  x  arises.  ^Once  x  is  determined  by  some  means,  o  can  be  obtained 
from  equation  (5)  and  m  from  solving  equation  (9).  Cohen  and  Woodward*  have 
published  tables  for  solution  of  the  equations,  and  we  have  prepared  a 
computer  program  in  BASIC  (see  the  Appendix)  for  numerical  solution. 

In  terms  of  the  lognormal  distribution  of  our  problem 


u' 

at 

in  d, 

u> 

0 

rm 

in  d  , 
o» 

u 

A 

m 

in  (d/dQ), 

A 

a 

a8 

m 

exp  (o),  and 

A 

dN 

m 

exp  (m) 

*  Cohen,  A.C.,  and  Woodward,  John.  Tables  of  Pearson-Lee-Fisher  Functions  of 
Singly  Truncated  Normal  Distribution.  Biometrics  9,  489-497  (1953). 
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A. 


ILLUSTRATIVE  EXAMPLE 


An  aerosol  sample  Injected  into  a  particle  counter  produced  the  data 
shown  in  table  1.  The  smallest  particle  resolved  was  0.15  um  in  diameter. 
One  can  see  that  the  data  are  truncated  just  below  the  mode  of  the 
distribution,  this  implies  that  neglect  of  the  truncation  could  produce  a 
significant  error.  In  fact,  the  estimates  obtained  by  the  blind  application 
of  routine  formulas  are  d„  «  0.200  and  o  *  1.20.  If  we  apply  the  Pearson-Lee 

algorithm  as  implemented  in  the  computer  program,  we  get  d^  =  0.186  and  o^  = 

1.26,  and  note  that  the  differences  in  the  estimates  are  in  the  expected 
directions. 

Table  1.  Experimental  Observations 


Point  ho.  Size  Range  Count 


1 

0.15 

- 

0.16 

375 

2 

0.16 

- 

0.17 

A08 

3 

0.17 

- 

0.18 

A28 

A 

0.18 

- 

0.19 

38A 

5 

0.19 

- 

0.20 

379 

6 

0.20 

- 

0.21 

333 

7 

0.21 

- 

0.22 

279 

8 

0,22 

- 

0.23 

207 

9 

0.23 

- 

0.2A 

203 

10 

0.2A 

- 

0.25 

139 

11 

0.25 

- 

0.26 

176 

12 

0.26 

- 

0.27 

113 

13 

0.27 

- 

0.28 

75 

1A 

0.28 

- 

0.29 

65 

15 

0.29 

_ 

0.30 

51 

5 .  SUMMARY 


Point  No. 

Size  Range 

Count 

16 

0.305  - 

0.330 

2  ’ 

17 

0.330  - 

0.355 

10 

18 

C.355  - 

0.380 

6 

19 

0.380  - 

0.A05 

3 

20 

0. AOS  - 

0.A30 

2 

21 

0.A30  - 

0.A55 

1 

22 

O 

-F~ 

L-n 

1 

0.A80 

0 

23 

0.A80  - 

0.505 

1 

2A 

0.505  - 

0.530 

0 

25 

0.530  - 

0.555 

0 

26 

0.555  - 

0.580 

0 

27 

0.580  - 

0.605 

0 

A  method  has  been  found  and  implemented  to  obtain  estimates  of  the 
parameters  of  truncated  normal  and  lognormal  distributions.  Emphasis  is 
placed  on  the  importance  of  recognizing  the  truncation  of  data  by  measuring 
devices  and  on  the  ease  with  which  the  problem  can  be  solved. 
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APPENDIX 


A  Computer  Program  for  Estimating  the  Parameters 
of  a  Truncated  Lognormal  Distribution 


A  computer  program  has  been  written  to  implement  the  Pearson-Lee 
method  of  estimating  the  parameters  of  a  normal  distribution  and  to  apply  the 
technique  to  data  from  a  lognormal  distribution.  The  program  is  written  in 
3ASIC  for  the  Tektronix  (TEK)  4051  Graphics  Calculator.  The  interactive  mode 
is  used  to  prompt  the  demand  terminal  user  for  the  proper  input  sequence. 
Because  the  TEK  4051  has  an  internal  tape  unit,  provision  is  made  for  data 
storage  and  retrieval  from  unit  33.  Statements  may  be  removed  from  unit  33 
and  used  with  other  terminals  or  host  computers.  The  program  also  calculates 
estimates  of  the  parameters  using  a  method  by  Cohen  which  uses  the  third 
moments  about  the  truncation  point.  Use  of  third  moments  is  believed  to 
introduce  inaccuracies  which  become  apparent  in  the  large  deviations  from  the 
Pearson-Lee  maximum  likelihood  estimates.  A  listing  of  the  program  appears  at 
the  end  of  the  appendix. 

Two  comments  are  needed  to  guide  the  program  user: 

1.  The  test  identification  can  be  a  string  of  72  characters. 

2.  The  size  requested  should  be  the  midpoint  of  the  size  range 
represented. 
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DISTRIBUTION  LIST  1 


Copies 


Names 


Copies 


Names 

CHEMICAL  SYSTEMS  LABORATORY 


ATTN:  DRDAR-CLB  1 

ATTN:  DRDAR-CLB -C  1 

ATTN:  DRDAR-CLB-PO  l 

ATTN:  DRDAR-CLB-PS  1 

ATTN:  DRDAR-CLB -R  1 

ATTN:  DRDAR-CLB-T  1 

ATTN:  DRDAR-CLB-TE  1 

ATTN:  DRDAR-CLC-B  1 

ATTN:  DRDAR-CLC-C  1 

ATTN:  DRDAR-CLC-E  1 

ATTN:  DRDAR-CLF  1 

ATTN:  DRDAR-CU-R  1 

ATTN:  DRDAR-CLJ-L  2 

ATTN:  DRDAR-CU-M  1 

ATTN:  DRDAR-CU-P  1 

ATTN:  DRDAR-CLN  1 

ATTN:  DRDAR-CLN-S  1 

ATTN:  DRDAR-CLN- ST  1 

ATTN:  DRDAR-CLN-TE  1 

ATTN:  DRDAR-CLT  1 

ATTN:  DRDAR-CLW-C  l 

ATTN:  DRDAR-CLW-P  1 

ATTN:  DRDAR-CLY-A  1 

AUTHORS  COPIES: 

ATTN:  DRDAR-CLY-A  11 

RECORD  COPY:  DRDAR-CLY-A  1 


DEPARTMENT  OF  DEFENSE 

Mr.  Joseph  F.  Campbell,  Jr. 

Office  of  Management  and  Budget /EXDP  1 

New  Executive  Office  Building  Rm  10002 
726  Jackson  Place,  N.W. 

Washington,  DC  20503 

Defense  Technical  Information  Center 
ATTN:  DTIC-DOA-2  12 

Caaeron  Station,  Building  5 
Alexandria,  VA  22314 


Comnander 
USASED,  USAINS00M 

ATTN:  IAFM-SED-III  l 

Fort  Meade,  MD  20755 

DEPARTMENT  OF  THE  ARMY 

HQDA 

ATTN:  DAMD-NCC  1 

WASH  DC  20310 

HQDA,  OCSA 

ATTN:  DACS-FM  (MG  R.  Anson)  1 

Room  1A871 ,  Pentagon 
Washington,  DC  20310 

Federal  Emergency  Management  Agency 
Office  of  Research/NPP 

ATTN:  David  W.  Benaen  1 

Washington,  DC  20472 

HQ  DA 

Office  of  the  Deputy  Chief  of  Staff  for 
Research,  Development  &  Acquisition 
ATTN:  DAMA-CSS-C  1 

Washington,  DO  20310 

Consaander 

HQ  7th  Medical  Command 

ATTN:  AEMPM  1 

APO  New  York  09403 

Consnander 
DARC0M,  STITEDR 

A3TN:  DRXST-STI  1 

Bex  48,  APO  New  York  09710 

Cagmander 

USASTCFEO 

ATTN:  MAJ  Mikeworth  1 

APO  San  Francisco  96328 
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Commander 

US  Army  Nuclear  4  Chemical  Agency 
ATTN:  MDNA-WE  1 

7500  Back lick  Rd,  Bldg  2073 
Springfield,  VA  22150 

Army  Research  Office 

ATTN:  DRXRO-CB  (Dr.  R.  Ghirardelli)  1 

P.0.  Box  12211 

Research  Triangle  Park,  NC  27709 
OFFICE  OF  THE  SURGEON  GENERAL 
Commander 

US  Army  Medical  Bioengineering  Research 
and  Development  Laboratory 
ATTN:  SGRD-UBD-AL ,  Bldg  568  1 

Fort  Detrick,  Frederick,  MD  21701 

Commander 

USA  Medical  Research  Institute  of 
Chemical  Defense 

ATTN:  SGRD-UV-L  1 

Aberdeen  Proving  Ground,  MD  21010 

US  ARMY  MATERIEL  DEVELOPMENT  AND 
READINESS  COMMAND 

Commander 
HQ,  DARCOM 

ATTN:  DRCED  (BG  Robinson)  1 

5001  Eisenhower  Ave 
Alexandria,  VA  22333 

Commander 

US  Array  Materiel  Development  and 
Readiness  Command 

ATTN:  DRCLDC  1 

ATTN:  DRCMT  1 

ATTN:  DRCSF-P  1 

5001  Eisenhower  Ave 
Alexandria,  VA  22333 


Coonander 

US  Array  Toxic  &  Hazardous  Materials  Agancy 


ATTN: 

DR3CTH-ES 

1 

ATTN: 

DRXTH-SE 

1 

ATTN: 

DRXTH-TE 

1 

Aberdeen  Proving  Ground,  MD  21010 
Commander 

US  Army  Foreign  Science  &  Technology  Center 
ATTN:  DRXST-MT3  1 

220  Seventh  St.,  NE 
Charlottesville,  VA  22901 

Director 

US  Array  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-MP  1 

ATTN:  DRXSY-CA  (Mr.  Metz)  1 

Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Array  Missile  Command 
Redstone  Scientific  Information  Center 
ATTN:  DRSMI-RPR  (Documents)  1 

Redstone  Arsenal,  AL  35809 

Director 

DARCCM  Field  Safety  Activity 

ATTN:  DRXOS-C  1 

Charlestown,  IN  47111 

Commander 

US  Army  Natick  Research  and  Development 


Laboratories 

ATTN:  DRDNA-UE  (AMEL)  1 
ATTN:  DRDNA-IC  1 
ATTN:  DRDNA-ICAA  1 
ATTN:  DRDNA-IM  1 
ATTN:  DRDNA-IB  1 
ATTN:  DRDNA-ITF  1 
Natick,  MA  01760 


Project  Manager  Smoke /Obscurants 

ATTN:  DRCPM-Sttt  1 

Aberdeen  Proving  Ground,  MD  21005 
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US  ARMY  ARMAMENT  RESEARCH  AMD 
DEVELOPMENT  COWAND 


Commander 

US  Army  Armament  Research  and 
Development  Command 

ATTN:  DRDAR-LCA-L  1 

ATTN:  DRQAR-LCE-C  1 

ATTN:  ORDAR-LCH-E  1 

ATTN:  DRDAR-LCM-SA  1 

ATTN:  DRDAR-LCU-CE  I 

ATTN:  DRDAR-NC  (COL  Lymn)  1 

ATTN:  DRDAR-SCA-A  1 

ATTN:  DPIWR-SCA-C  1 

ATTN:  DRDAR-SCA-T  1 

ATTN:  DRDAR-SCM  1 

ATTN:  DRDAR-SCS  1 

ATTN:  DRDAR-TDC  (Dr.  D.  Gyorog)  1 

ATTN:  DRDAR-TSS  2 

ATTN:  DRCPM-CAWS-AM  1 

ATTN:  DRCPM-CAHS-SI  1 

Dover,  NJ  07S01 


Commander 

ARRADOOM 

ATTN:  DRDAR-QAC-E  1 

Aberdeen  Proving  Ground,  MD  21010 

Director 

USA  Ballistic  Research  Laboratory 
ARRADOOM 

ATTN:  DRDAR-TS8-S  1 

Aberdeeu  Proving  Ground,  MD  21005 

US  Army  Armament  Research  and 
Development  Conraand 
Resident  Operations  Office 
ATTN:  DRDAR-TSE-OA  (Robert  Thresher)  1 

National  Space  Technology  Laboratories 
NSTL  Station,  MS  39529 


US  ARMY  ARMAMENT  MATERIEL  READINESS 
COfiAND 

Commander 

US  Army  Armament  Materiel  Readiness  Ccranand 


ATTN:  DRSAR-ASN  1 

ATTN:  DRSAR-IRW  1 

ATTN:  DRSAR-ISE  1 

ATTN:  DRSAR-PE  1 


Rock  Island,  IL  61299 

Ccomander 
USA  ARRCOM 

ATTN:  DRSAR-MAY-M  1 

ATTN:  DRSAR-MAS-C  1 

ATTN:  SARTE  l 

Aberdeen  Proving  Ground,  MD  21010 

Coemander 

US  A. my  Dugway  Proving  Ground 

ATTN:  Technical  Library  (Docu  Sect)  l 

Dugway,  UT  84022 

US  ARMY  TRAINING  &  DOCTRINE  COMMAND 
Ccnmandant 

US  Army  Infantry  School 

ATTN:  CTDD,  CSD,  NBC  Branch  1 

Fort  Betming,  GA  31905 

Conmandant 

US  Army  Missile  &  Munitions  Center 
and  School 

ATTN:  AXSK-TME  1 

Redstone  Arsenal,  AL  35809 

Commander 

US  Army  Logistics  Center 

ATTN:  ATCL-*C  1 

Fort  Lee,  VA  23801 
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Conmand  ant 

US  Array  Chemical  School 

ATTN:  ATZN-CM-C  l 

ATTN:  ATZN-CM-AD  2 

ATTN:  ATZN-CM-TPC  2 

Fort  McClellan,  AL  36205 

Conmander 

USAAVNC 

ATTN:  AIZQ-D-MS  1 

Fort  Rucker,  AL  36362 

Commander 

US  Array  Infantry  Center 

ATTN:  ATSH-CD-MS-C  1 

Fort  Benning,  GA  31905 

Conmander 

US  Array  Infantry  Center 
Directorate  of  Plana  &  Training 
ATTN:  ATZB-DPT-PO-NBC  1 

Fort  Benning,  GA  31905 

Conmander 

USA  Training  and  Doctrine  Conmand 

ATTN:  ATCD-N  1 

Fort  Monroe,  VA  23651 

Conmander 

US  Array  Armor  Center 

ATTN:  ATZK-CD-MS  1 

ATTN:  ATZK-DPT-PO-C  1 

Fort  Knox,  KY  40121 

US  ARMY  TEST  S  EVALUATION  COMMAND 

Conmander 

US  Array  Test  &  Evaluation  Comaand 
ATTN:  DRSTE-CM-F 

Aberdeen  Proving  Ground,  H)  21005 

DEPARTMENT  OF  TV-  NAVY 

Chief,  Bureau  of  Medicine  &  Surgery 
Department  of  the  Navy 
ATTN:  MED  3C33 

Washington,  DC  20372 


Project  Manager 

Theatre  Nuclear  Warfare  Project  Office 
ATTN:  TN-09C  1 

Navy  Department 
Washington,  DC  20360 

Chief  of  Naval  Research 
ATTN:  Code  441  1 

800  N.  Quincy  Street 
Arlington,  VA  22217 

Conmander 

Naval  Facilities  Conmand 
ATTN:  FAC-032  1 

ATTN:  FAC-062  1 

200  Stovall  St 
Alexandria,  VA  22332 

Conmander 

Naval  Surface  Weapons  Center 
Code  G51  1 

Dahlgren,  VA  22443 

Commanding  Officer 
Naval  Weapons  Support  Center 
Applied  Sciences  Department 
ATTN:  Code  50C,  Bldg  190  1 

Crane,  IN  47522 

Conmander 

Naval  Air  Development  Center 
ATTN:  Code  2012  (Dr.  Robert  Helnbold)  1 

Warminster,  PA  18974 

US  MARINE  CORPS 

Conmand ing  General 
Marine  Corps  Development  and 
1  Education  Conmand 

ATTN:  Fire  Power  Division,  D091  1 

(^jantico,  VA  22134 

DEPARTMENT  OF  THE  AIR  FORCE 

Department  of  the  Air  Force 
l  Headquarters  Foreign  Technology  Division 

ATTN:  TQTR  1 

Wright-Patterson  AFB,  CXI  45433 
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HQ  FTD/NIIS  1 

Wright-Patterson  AFB,  OH  45433 

ASD/AESD  1 

Wright-Patterson  AFB,  OH  45433 

AFAMRL/HE 

ATTN:  Dr.  Clyde  Reploggle  1 

Wright-Patterson  AFB,  OH  45433 

AFAMRL/TS 

ATTN:  COL  Johnson  1 

Wright-Patterson  AFB,  OH  45433 

HQ  AFSC/S02  1 

ATTN:  CPT  D.  Riediger 

.Andrews  AFB,  MD  20334 

HQ,  AFSC/SDNE  1 

Andrews  AFB,  MD  20334 

HQ,  AFSC/SGB  1 

Andrews  AFB,  DC  20334 

HQ,  NORAD 

ATTN:  J-3TU  1 

Peterson  AFB,  CO  30914 

USAF  TAWC/THL  1 

Eglin  AFB,  FL  32542 

AFATL/DLV  1 

Eglin  AFB,  FL  32542 

USAF  SC 

ATTN:  AD/YQ  1 

ATTN:  AD/YQO  (MAJ  Owens)  1 

Eglin  AFB,  FL  32542 

AD/XRO  1 

Eglin  AFB,  FL  32542 


USAFSAM/VN 

Deputy  for  Chemical  Defense 

ATTN:  Dr.  F.  Wesley  Baungardner  l 

Brooks  AFB,  TX  78235 


AH3/RDTK 

ATTN:  LTC  T.  Kingery  1 

Brooks  AFB,  TX  78235 

OUTSIDE  AGENCIES 

Battelle,  Columbus  Laboratories 
ATTN:  TACTEC  1 

505  King  Avenue 
Columbus,  OH  43201 

Toxicology  Information  Center,  JH  652 
National  Research  Council  1 

2101  Constitution  Ave.,  NW 
Washington,  DC  20418 

US  Public  Health  Service 

Center  for  Disease  Control 

ATTN:  Lewis  Webb,  Jr.  1 

Building  4,  Room  232 

Atlanta,  GA  30333 

Director 

Central  Intelligence  Agency 

ATTN:  AMR/ORO/DD/S&T  1 

Washington,  DC  20505 

ADDITIONAL  ADDRESSEES 

Commandant 

Academy  of  Health  Sciences,  US  Army 
ATTN:  HSHA-CDH  1 

ATTN:  HSHA-IPM  2 

Fort  Sam  Houston,  TX  78234 

Commander 

217th  Chemical  Detachment 

ATTN:  AFVL-CD  1 

Fort  Knox,  KY  40121 

Headquarters 

US  Army  Medical  Research  and 
Development  Catmand 

ATTN:  SGRD-RMS  1 

Fort  Detrick,  MD  21701 
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Ccmnander 

US  Army  Envirotmental  Hygiene  Agency 
ATTN:  Librarian,  Bldg  2100 
Aberdeen  Proving  Ground,  MD  21010 

Stimaon  Library  (Documents)  1 

Academy  of  Health  Sciences,  US  Army 
Bldg.  2 8 AO 

Fort  Sam  Houston,  TX  78234 
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